Influence of fluctuations in actin structure on myosin V step size 
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We study the influence of disorder in the helical actin structure on the myosin V step size, 
predicted from the elastic lever arm model. We show that fluctuations of ±5 degrees per actin 
subunit, as proposed by Egelman et al., significantly alter the distribution of step sizes and improve 
the agreement with experimental data. 



Introduction 

Myosin V is a motor protein from the myosin superfam- 
ily, involved in various intracellular transport processes 
[l| . It is a processive motor |2j , which means that a sin- 
gle protein molecule can transport cargoes along actin 
filaments over distances of several micrometers. It is 
dimeric, consisting of two identical heads, each attached 
to a lever arm, joined with each other through a tail do- 
main. It achieves processivity by alternately binding its 
two heads to an actin filament and thus walking in a 
hand-over-hand fashion Q. Its step size is roughly de- 
termined by the periodicity of the actin filament, which 
is about 36 ran. 

In comparison with muscle myosin (myosin II), two 
major adaptations are found in myosin V: a longer lever- 
arm, measuring around 26 nm |^ and a slower release of 
ADP. The duty cycle of each head is otherwise similar 
0: the head binds to actin in the ADP.Pi state, first 
releases Pi and performs a large conformational change 
(power stroke), then releases ADP and performs a smaller 
conformational change, and finally binds a new ATP 
molecule and detaches from actin. 

In a recent article @, we have developed a model for 
myosin V, based on the elasticity of the lever arm. It 
assumes that the lever arm is stiffly anchored in the 
head, but the direction depends on the chemical state 
of that head. The distal ends of the lever arms are con- 
nected together through a flexible hinge, which repre- 
sents the only means of "communication" between the 
heads. By calculating the bending energies in all rele- 
vant dimer states and its influence on transition rates, 
we have shown that the elastic lever-arm model explains 
the coordinated hand-over-hand motility by showing that 
the lead head is not likely to bind to actin before the 
trail head undergoes the major power stroke and that 
it cannot commit its power stroke before the trail head 
unbinds. The model also quantitatively reproduces the 
measured force- velocity relations and shows how the run 
length (the average distance a motor runs processively 
before it dissociates from actin) could be used to deter- 
mine some kinetic rates. Although the original model 
reproduces the step size corresponding to about one half- 
turn of the actin helix (13 actin subunits), there is a small 
but significant deviation between the predicted distribu- 
tion of step sizes (mainly on the 13th subunit, with a side 
peak on the 11th) and the statistics obtained from elec- 
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FIG. 1: The regular 13/6 helical model of an actin filament 
(top), and the 28/13 helix with ±5° fluctuations per subunit 
(bottom). The subunits are shown blurred, with a relative 
intensity corresponding to their statistical weight. The orien- 
tation of the central subunit (denoted as "0") is always kept 
fixed. 



tron microscopy (EM) studies Q, which show the main 
peak on the 13th subunit and two significant side peaks 
on the 11th and the 15th subunit, whereby the 15th is 
stronger than the 11th. 

The purpose of the current paper is to extend the 
present model by including a more precise description 
of the actin structure, which could potentially explain 
the current deviation. The central modification will be 
to take into account torsional fluctuations in the actin 
structure. So far we have used the assumption that the 
actin monomers are arranged on a helix with a periodic- 
ity of 13 actin subunits (also called a 13/6 helix, because 
the rotation of each next subunit is ^ x 360°. However, 
it has been known for some time that the actin filaments 
have a variable, fluctuating twist The nature and 
the dynamics of these fluctuations are not yet entirely 
understood, but it seems that the azimuthal orientation 
of each subunit can fluctuate about ±5° with respect to 
its neighbors 0- There is some evidence that the dy- 
namics of these fluctuations is rather slow 0, |n| , with 
a characteristic time of the order of seconds. This would 
imply that the twist has discrete states with a consider- 
able energy barrier between them. The average twist per 
subunit was found to be about 167°, slightly more than 
in the 13/6 helix, where it would be 166.15° (Fig. HJ. 

Based on many experimental and theoretical studies 
IE 13 113 ■ the following picture of the working cy- 
cle of dimeric myosin V has now found broad consensus 
(Fig. EJ. The major part of the dimeric cycle is spent in 
the state with both heads binding ADP. The trail head 
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TABLE I: Model parameters 



Lever arm length 


L 


26 nm for 6IQ (WT) 


Lever arm start (radial) 


R 


8 nm 


Lever arm start (longitudinal) 


^ADP.Pi 


nm 




<5adp 


3.5 nm 


Angle ADP.Pi 


</>ADP.Pi 


115° 


Angle ADP 


0ADP 


50° 


Bending modulus 


EI 


1500 pN nm 2 


Average helix twist per subunit 


0° 


167.14° 


Twist fluctuations 


0' 


5° 


Helix rise per subunit 


a 


2.75 nm 



is in the post-powerstroke and the lead head likely in the 
pre-powerstroke state 0, 0, \^ > although the latter has 
still been controversial recently (cf. ref. ^J). Then the 
trail head releases ADP, binds a new ATP molecule and 
unbinds from actin. With the trail head released, the 
lead head is now free to undergo a power-stroke, where- 
upon the former trail head binds in the lead position 
(one step ahead in the direction of motion) and releases 
Pi. This brings the dimer to the original state, however 
one step further toward the actin plus end and with one 
ATP molecule hydrolyzed. 



Results 

To calculate the step size, which is the main objective 
of this article, the relevant transition is the binding of the 
lead head (in the ADP.Pi state) while the trail head is 
in the forward-leaning ADP conformation. Let us denote 
the binding rate to the site i (counted from the trail head, 
which is bound to site 0) with k + A(i) and the rate of the 
reverse (unbinding) process as k^AH)- From the detailed 
balance we know that the ratio of these rates has to obey 



e +A(i) 



K +A 
1,0 



exp 



AU + FAx 



(1) 



where AU denotes the increase in elastic energy upon 
binding and FAx denotes the work done against the 
applied load. fc° A and k°_ A represent the binding and 
unbinding rates of a head that binds without any elas- 
tic distortions in the lever arm. Because the persistence 
length of the lever arms is significantly longer than their 
actual length, we neglect thermal fluctuations and the 
corresponding entropic contributions to the total free en- 
ergy. In the following we will also concentrate on the case 
of an unloaded motor (F = 0). In this case the bending 
energy in the initial state is and AU is equal to the 
elastic energy in the double-bound state. All model pa- 
rameters are summarized in Table [I] and Fig. [21 The 
geometric parameters were obtained from published EM 
studies 0,0. The bending modulus of the lever arm (EI) 
was estimated based on a lower boundary provided by 



the fact that a myosin V molecule is able toperform reg- 
ular steps against loads of at least 1.8 pN |16| . 

The elastic energy of a certain dimer conformation is 
calculated in the following way (which is described in 
more detail in Ref.|(|). A lever arm starting point is 
positioned at 



(2) 




and its initial tangent is 

^ cos(0) 

sin(0) sin(6*j) 
^ sin(</>) cos(6'i) 



(3) 



Other conditions state that the length of the lever arm 
from the starting point to the joint has to be L and that 
the endpoints of both lever arms coincide. The bending 
energy is then given as 



U 



2 

E 

3=1 



EI 



ds 



ds 



(4) 



where EI denotes the bending modulus and 4^ the local 
curvature of a lever arm. The summation index j runs 
over both lever arms. The shape of the dimeric molecule 
is determined numerically in a way that minimizes U. 

There are two possibilities concerning the reversibility 
of the binding process. One possibility is that the lead 
head binds to a site and does not detach before the next 
power stroke. In this case, the probability that it binds 
to the site i is 



(5) 



If, however, the lead head in the ADP.Pi (weakly bound) 
state is allowed to detach and re-attach several times, 
the probability to find it on the site i is given by the 
equilibrium distribution, 



Pi = 



k+A(i)/k-A(i) 

k +AU)/k- A ( 3 ) 



(6) 



So far, no experimental data are available about the 
strain dependence of the detachment rate, k_A- How- 
ever, as with most binding processes one can expect that 
the activation point for the binding transition is closer to 
the bound state, and therefore the strain-dependence of 
k-A is weaker than that of k + A- We therefore neglect the 
strain-dependence of detachment rates and assume that 
the Boltzmann factor in Eq. only influences k + A(i)- 
If the binding process is reversible, Pi (Eq. [BJ) becomes 
independent of this assumption in any case. 

To describe the fluctuations in the actin helix, we will 
use the coefficients a,, which can assume the values -1, 
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Direction of motion 



FIG. 2: The most likely duty cycle of the dimeric myosin V. A head in the ADP state undergoes the power-stroke, then the 
second head containing ADP.Pi binds in the lead position, it releases Pi, then the rear head releases ADP (time limiting), binds 
a new ATP molecule and detaches from actin and the cycle repeats with exchanged roles of the two heads, however with the 
dimer having proceeded about one actin period further. 




FIG. 3: The geometry of a dimeric myosin V molecule, after 
the lead head has bound to a site 13 subunits in front of the 
trail head. 

or 1. This corresponds to the assumption that each 
subunit can have three orientations relative to its left 
neighbor [To| . The twist between the subunits i — 1 and i 
is then 0j — = 9° + ai9' , and the total twist of subunit 

1 relative to subunit is 9i — i9° + (ai + a 2 + . . . a,)0'. 

Another important question is the dynamics of the 
actin fluctuations. Let us denote the probability that 
the helix is in the state with angles 9° + a\9' , 9° + a^9' ', 
. . . with 



ai,a2,..., 




If the fluctuations are fast in comparison with the at- 
tachment rate, the probability that the head binds to site 
i is proportional to the attachment rate in each configu- 
ration, weighted by the probability of that configuration 
and summed over all possible configurations: 

V P p -Ut{0°i+0'(a 1 +a 2 + ...a z ))/k B T 

/ , 1 Oi,02,---,Oi c 
p/ 01,02, •••«» 

i -£ £ p aia2 ^ aie -U ] (0° 3 +e'(a 1+ a 2 + ...a J ))/k B T 

(8) 

Here Ui (9) denotes the elastic energy in a state where the 
lead head is bound i subunits in front of the trail head 
and the helical twist between these two subunits is 9. A 
simplified expression for Pf is derived in the Appendix. 



If, on the other hand, the fluctuations are slow, the 
probability is given by the ensemble-average of all helix 
conformations: 

p -U,(e°i+e'(a 1 +a 2 + ...a t ))/k B T 

p? = V p — 

r t ^oi,o 2 ,...,o„y, e -U J (e a j+6>(a 1 +a 2 + ...a J ))/k B T 

01,02,... . 

3 

(9) 

Here n denotes the maximum index of a site that still has 
a non-negligible binding probability. In the following, we 
will calculate the results for both scenarios, even though 
we consider the slowly fluctuating scenario more realistic. 

A numerical evaluation of the probability distributions 
P" and P( for four different lever arm lengths, corre- 
sponding to mutants with 2, 4, 6 (wild type) and 8 IQ 
motifs in the lever arm shows a notable difference in com- 
parison with results on the 13/6 helix (Fig.QJ. 

To check the effects of uncertainties in other model 
parameters on our results, we have tried three modifica- 
tions: we have introduced an additional azimuthal com- 
ponent of the power-stroke (meaning that after Pi release, 
the lever arm makes an additional swing of 20° to the 
right, while its starting point remains the same), we have 
changed the lever arm angle <fi in the pre-powerstroke 
state and we have changed the lever arm stiffness. The re- 
sults are summarized in TablelTTl None of the tested mod- 
ifications makes an improvement with regard to agree- 
ment with experimental data as compared to the basic 
model. 

Another assumption we made in our model is that the 
lever arm anchoring is stiff in every state. While this 
assumption is justified for long lever arms, the effect of 
compliance in the converter domain could play an impor- 
tant role for short lever arms. Results on skeletal muscle 
myosin (with short lever arms) indeed show that at least 
half of the compliance is in the converter domain [l7| . 
Therefore some care has to be taken when interpreting 
the results for 2IQ lever arms. 
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TABLE II: Average step size for the basic model and the following variations: (i) additional azimuthal power stroke, (ii) 
different lever arm angle in the pre-powerstroke conformation (<f>ADP.Pi = 135° instead of 115°) and (iii) a softer lever-arm 
(EI — 500pNnm 2 instead of 1500 pNnm 2 ). All data are based on the slow-fluctuation scenario. 



Lever 


arm length basic model A9 = 


20° (f>ADP.Pi 


= 135° EI = 500 pN ran 2 


2IQ 


2.0 


2.28 


2.6 


2.1 


4IQ 


8.1 


4.5 


11.7 


7.9 


6IQ 


13.4 


14.1 


14.0 


13.9 


8IQ 


14.1 


15.1 


14.7 


13.2 



(A) 



(B) 
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6IQ 
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i i i i i i i i 
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8IQ 
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Step size (actin subunits) 
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4IQ 



i i i , 



6IQ 



i i i i i i i 



8IQ 







I I 
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Step size (actin subunits) 



Discussion 

The results for the slowly fluctuating scenario show a 
significantly improved agreement with available experi- 
ments than the previous model, based on a stiff 13/6 
helix. In particular, the observation that the probabil- 
ity for 15 subunit steps is about twice as high as for 11 
subunit steps is well reproduced. The average step size 
of 13.4 actin subunits (on a helix with a half-pitch of 
14) also means that the molecule makes a left-handed 
rotation of about 8° per step, or 0.2° per nm traveled. 
This means that a freely walking myosin V makes one 
revolution around the actin filament every 1.7 /im, well 
consistent with a value of 2.2 /im, measured by Ali and 
co-workers [l8j . Note that their interpretation (based on 
the assumption of a 13/6 helix) is that the steps repre- 
sent a mixture of 11 and 13 subunit lengths. However, 
on a 28/13 helix, data become consistent with a mixture 
including 15 subunit steps, which are even more frequent 
than 11 subunit steps. 

The third experiment our model should be tested 
against is the dependence of the step size on the lever arm 
length. Purcell and co-workers 0] used optical tweezers 
to determine the step size of different dimeric constructs. 
For wild- type myosin V (6 IQ domains), they measured 
a step size of 35nm. For 4IQ mutants they measured 
24nm or 8.7 actin subunits, and for HQ mutants 5nm 
or 2 subunits. Also in this respect the current model 
(predicting 13.4, 8 and 2 subunit steps, respectively) no- 
ticeably improves the agreement with experimental data. 
The previous model, based on a stiff helix, predicted a 
more abrupt step size change between 2 and 4IQ mutants. 



Appendix 



FIG. 4: The probabilities that the lead head binds to site i 
in the model with slow fluctuations (P/, A) and fast fluctu- 
ations (P/ , B). The four diagrams are for different lever arm 
lengths (2IQ: lOnm; 4IQ: 18nm; 6IQ: 26nm; 8IQ: 34nm). The 
vertical lines in each diagram show the average position (i). 
The dashed lines represent results on a rigid 13/6 helix, as 
used in previous calculations. 



In Eq. © the summation over coefficients di, . . . , a% 
can be carried out while keeping their sum constant (si). 
This way we obtain the probability that the subunit i is 
oriented at angle 9°i + 9' Si\ 



P 



Pn 



i . , ^y^j '■ < : ) . ii ■> < i ... 

a± ,a,2 ,. . .,fli;ai +0-2 + - --+0-1— Si 



k;\si\<k<(i+\si\)/2 




(10) 



5 



Then expression JSJ can be simplified to Acknowledgment 

E Pi s , e -Ui{e°i+e' St )/k B T 

P f = . (ii) 

E E Pj, Sj e- u ^ e ^+ e '^y k B T 

For large values of i the central limit theorem can be This work was supported by the Slovenian Office of 

applied and the distribution P^ Si becomes Gaussian. Science (Grants No. Zl-4509-0106-02 and PO-0524-0106). 
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